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1. Introduction 

Study of heavy quarkona is important to understand the quark gluon plasma (QGP) formation 
in heavy ion collision experiments e.g. the RHIC experiment at Brookhaven National Laboratory. 
Recent studies on the spectral function of charmonium above Tc suggest that hadronic excitations 
corresponding to J/y may survive in the deconfinement phase till relatively high temperature 
^ |3[ m5[, ^]. Such results of strongly interacting QGP may affect the scenario of J/y suppression 
[0, ^]. Therefore the determination of the accurate dissociation temperature of 7/i/Ais required for 
many phenomenological studies. 

In order to study temporal correlators at finite temperature, one has to take account of a charac- 
teristic contribution caused by the finite temporal extent. Usually a meson correlator is interpreted 
by a diagram with quark and anti-quark propagators like it is sketched in Fig. |l|(a). However in the 
case of a system with a finite temporal extent, a wraparound (scattering) contribution, as shown in 
Fig. [l](b), has to be included as well. The effects are also known even at zero temperature case e.g. 



the study of two pion correlators [90, and pentaquark correlator [10|. In the case, the contribution 



is usually removed using e.g. Dirichlet boundary conditions ||10|] or appropriate analyses [Q]. On 
the other hand, at finite temperature, the latter contribution is genuine, and it provides a constant 
contribution to the correlator as a zero energy mode. Due to the contribution, correlators of such 
meson-like operators, e.g. charmonia, may be drastically changed just after the deconfinement 
transition. 

In this proceeding, we show how the constant mode affects the thermal effect in the charmo- 
nium correlators or the spectral functions. Details of the study has already been reported in the 



published paper [11]. 



2. Constant mode in free quark case 

2.1 Meson-like correlators in the free quark case 

In order to investigate the constant contribution, we first consider the free quark case of QCD, 
in which the constant contribution can be easily calculated. Here we define the meson-like corre- 
lators with quark bilinear operators, Or{x,t) = q{x,t)Tq{x,t), in the free quark case this gives for 
the correlator, 

C(0 = L(Cr(x,Oo[.(0,0)), (2.1) 

X 

where T are appropriate 4x4 matrices, i.e. 75, 7, 1, and 775 for pseudoscalar(Ps), vector(V), 
scalar(S), and axialvector(Av) channels respectively. In this paper we calculate correlators for 
degenerate quark masses, constructed from bilinear operators, with vanishing spatial momentum 
as the simplest case, which is aiming at studies of charmonium at finite temperature discussed latter. 
The spectral function of the correlator is defined by 



poo 

C{t) = / d(opr{(o)K{(o,t), 

JO 



cosh(a)(^-0) 

Ki(o,t) = / ^^ (2.2) 

sinh {(Of) 



and has been calculated in Ref. [12, 13 1 in the high temperature Umit. 
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Figure 1: A sketch of quark line diagrams for a meson-like correlator in a system with finite temporal extent. 
Vertical lines show the boundaries in temporal direction. 



2.2 Free quark calculations on a lattice 



Here we present numerical calculations of the meson-like correlators for free quarks on a 
lattice. The calculations are performed on an isotropic x A/^ = 16^ x 32 lattice. The free quark 
is described by the Wilson quark action with a bare quark mass of niq = 0.2. 

Figure ^ (left panel) shows the effective masses of the correlators defined above. The effective 
mass is defined by correlator at successive time slices. 



C{t) 



cosh ['Meff(0(f -0] 



C{t + \) cosh[mgff(? 



1)] 



(2.3) 



The effective masses of the meson-like correlators with free quarks should approach the energy of 
the two quark state without momentum for S-wave states (Ps and V channels) and with a minimum 
momentum for P-wave states (S and Av channels) except for the zero energy mode. In the Ps and V 
channels, the zero relative momentum processes have no constant contribution and their effective 
masses approach the expected values. In the S and Av channels, on the other hand, their effective 
masses approach the zero energy level rather than the energy of two quark states. 

2.3 An analysis to avoid the constant contribution 

Next we consider how to estimate the effect of the constant mode. Of course the cosh -i- 
constant fit enables us to do so, but here we present the method without a fit analysis. The method 
uses the midpoint subtracted correlator. 



C{t)=C{t)-C{Nt/2). 
One can also define the effective mass of the subtracted correlator. 



(2.4) 



C(0 sinh' 






C{t + l) sinh^ 


KffW(T-^-i); 



(2.5) 



Figure Q (right panel) shows the effective mass m^^(f), defined in Eq. (2.5), from the free 
quark results for Ns = 16. The effective masses are equivalent to the usual effective mass shown 
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Figure 2: (Left) Effective masses of meson-like correlators for each channel. These are free quark calcula- 
tions on a 16^ X 32 lattice. 



in Fig. ^ (left panel) except for the effects of the constant mode. The expected energies of the 
lowest two quark states for the S-wave and the P-wave states are shown in both figures of Fig. |2l In 
contrast to the case of the usual effective masses, the subtracted effective masses m^^{t) approach 
the expected values even in the P-wave states. The analysis to avoid the constant contribution works 
well at least in the free quark case. 

The method using the midpoint subtracted correlators can also be applied to studies of spectral 
functions Pr(w) by a modification of the kernel, 

C{t) = [ d(opr{(o)K^^^{(o,t), (2.6) 







,,,u, , 2sinh2(f (^-0) 

= . • (2-7) 

smh (ft>Y) 

By using the alternative kernel K^^^{(0,t), it is possible to extract the spectral function without the 
contribution in ft) ^ T from e.g. the Maximum Entropy Method. 

3. Quenched QCD case at finite temperature 
3.1 Lattice setup 

In this section we demonstrate effect of the constant mode and the analysis to avoid it in 
quenched QCD at finite temperature. Here we calculate the charmonium correlators, which are 



defined like in Sect. ^ at the charm quark mass. 

The gauge configurations have been generated by an standard plaquette gauge action with 
a lattice gauge coupling constant, j3 = 6.10 and a bare anisotropy parameter jc = 3.2108. The 



definition of the action and parameters are the same as that adopted in Ref. [|14j]. Lattice sizes 
are 20^ x A^, where A^j = 160 at T = and A^, = 32, 26 and 20 at T > 0. The lattice spacings are 
l/a^ = 2.030(13) GeV and l/ot = 8.12(5) GeV which are determined from the hadronic radius 
ro = 0.5 fm. The physical volume size is about (2fm)^ and temperatures at Nt = 32, 26 and 20 ai^e 
T = OMTc, imTc and lAOTc respectively. 
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Figure 3: (Left) Effective masses from usual charmonium correlators. (Right) Effective masses from the 
midpoint subtracted correlators. These are results with local operators for each channel in quenched QCD at 
finite temperature. 



For the quark fields, we adopt an 0{a) improved Wilson quark action with tadpole improved 
tree level clover coefficients. Although the definition of the quark action is the same as in Ref. [p^], 
we adopt a different choice of the Wilson parameter r = 1 to suppress effects of lattice artifacts in 
higher states of charmonium [12]. 



3.2 Finite temperature results 

First we show the usual effective masses, m^f^{t) at finite temperature in Fig. ^(Left panel). 
Since these are results obtained with local operators, there is no plateau region at any temperature 
or channel we chose. It can be expected from zero temperature calculation and the effective masses 
show reasonable behavior below Tc, at which the charmonium correlators does not show large 
difference from the zero temperature results [§, ^, ^]. Just above T^, we see drastic changes in 
the P-wave state channels. The changes may cause the dissolution of Xc states just above Tc as 
already reported in Ref. [|3|, |5|, |6|]. On the other hand, the S-wave states, i.e. the Ps and V channels, 
show small changes up to lATc. This temperature dependences are shown in Fig. ^ (left panel) 
for each channel. When we compare the effective masses at the same time slices, we find no large 
change up to 1 ATc in both channels as well. The results are consistent with previous lattice studies 
of charmonium spectral functions [0, ^ ^, ^, ^, in which the spectral functions for the Ps and V 
channels in the deconfinement phase (but at not so high temperature) show a peak structure around 
the T]c and masses like the result at zero temperature. 



3.3 Midpoint subtracted correlator analysis 

Figure ^ (right panel) shows the results of m^^{t) in quenched QCD at finite temperature. In 
the P-wave states, the drastic changes of the usual effective masses wig^f (?) are absent in the effec- 
tive masses from the midpoint subtracted correlators m^^(f). Furthermore the similar behaviors 
of the effective masses hold till \ATc as in the case of the S-wave states. The results indicate that 
the drastic changes of the usual correlators (and effective masses) are caused only by the constant 
mode to the correlators. The situation is very similar to the free quark case discussed in Sect. ^ 
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Figure 4: (Left) Effective masses from usual charmonium correlators. (Right) Effective masses from the 
midpoint subtracted correlators. These are results with the spatially extended operators for each channel in 
quenched QCD at finite temperature. 



3.4 Results with spatially extended operators 

Finally we present results of meson correlators with spatially extended operators. The spatially 
extended operators are constructed with a smearing function which yields a spatial distribution 
of quark (and anti-quark) source(s). The spatially extended operators are defined by OY{x,t) = 
Y,y^{y)q{x — y,t)Yq{x,t) with a smearing function 0(x) in Coulomb gauge. In this calculation the 
spatially extended operators are adopted only in the source operators, the sink operators are local in 
any cases. The smearing function is the same as that in Ref. [Q], i.e. 0(x) = exp(— A[x|^) where A 
and P are parameters determined by a matching with the charmonium wave function as A = 0.2275 
and P = 1.258. Figure ^ (left panel) is the result for the usual effective masses with the spatially 
extended operators. Other conditions are the same as in Fig. ^. Below Tc, in contrast to the case 
of local operators the effective masses reach a plateau due to the larger overlap with the lowest 
state. Above T^- the effective masses of the P-wave states change more than in the case of local 
operators. Because the constant contribution is enhanced by the smearing function more than the 
other contributions in the case. Figure ^ (right panel) is the same figure as Fig. ^, but with spatially 
extended operators. We find similar behavior for both cases of operators. 

4. Conclusion 

Let us discuss details on Xc states using the results of the previous sections. In this study 
we find that the drastic changes of the (usual) correlators in the P-wave states are caused by the 
constant contribution, while the correlators without the constant mode yield small change till, at 
least, lATc even in the P-wave states. The changes are of similar size as that of the S-wave states 
as one can see in Fig. |3| and Fig. 

Although the constant mode coming from the scattering process is important to investigate 
some transport phenomena, it does not affect the spectral function in the (O ^ T region for the 
charmonium correlators with zero spatial momentum. Therefore the dissolution of Xc states just 
above Tc as discussed in e.g. Ref. |j3|, |5[, ^ might be misleading. When there is the constant contri- 
bution, the hadronic state is no longer the lowest state. The situation provides some difficulties in 
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the analysis to reconstruct the spectral function in the higher energy part (O^T even if the MEM 



is applied [|15|]. 

In principle the MEM analysis can extract the correct spectral function even if the constant 
contribution exists, however it is difficult to reproduce the non-lowest part of the spectral function 
correctly by the MEM analysis at finite temperature with typical statistics. In order to investigate 
the properties of Xc states above r^., a careful analysis taking the constant contribution into account 
is necessary. A correct MEM analysis of the correlator with and without the constant contribution 
should give the same results, up to the low energy region (o <^T. Such a comparison can be easily 



performed by the MEM analysis with the alternative kernel of Eq. (2.7) for midpoint subtracted 
correlators. Failure to do so signals unreliable MEM results. 
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